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ABSTRACT 

The so-called "acoustic revival mechanism" of core-collapse supernova pro- 
posed recently by the Arizona group is an interesting new possibility. Aiming 
to understand the elementary processes involved in the mechanism, we have cal- 
culated the eigen frequencies and eigen functions for the g-mode oscillations of 
a non-rotating proto neutron star. The g-modes have relatively small eigen fre- 
quencies. Since there is a convective region inside the proto neutron star, the 
eigen functions are confined either to the exterior or interior of this evanescent 
region although some fraction of eigen functions penetrates it because of their 
long wave lengths. The possible excitation of these modes by the standing ac- 
cretion shock instability, or SASI, is discussed based on these eigen functions. 
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We have formulated the forced oscillations of g-modes by the external pressure 
perturbations exerted on the proto neutron star surface. The driving pressure 
fluctuations have been adopted from our previous computations of the axisym- 
metric SASI in the non-linear regime. We have paid particular attention to low 
I modes, since these are the modes that are dominant in SASI and that the Ari- 
zona group claimed played an important role in their acoustic revival scenario. 
Here i is the index of the spherical harmonic functions, Y™ . Although the fre- 
quency spectrum of the non-linear SASI is broadened substantially by non-linear 
couplings, the typical frequency is still much smaller than those of g-modes, the 
fact leading to a severe impedance mismatch. As a result, the excitations of 
various g-modes are rather inefficient and the energy of the saturated g-modes is 
~ 10 50 erg or smaller, with the #2 mode being the largest in our model. Here the 
#2 mode has two radial nodes and is confined to the interior of the convection 
region. The energy transfer rate from the g-modes to out-going sound waves is 
estimated from the growth of the g-modes and found to be ~ 10 51 erg/s in the 
model studied in this paper. 

Subject headings: supernovae: collapse — neutrinos — hydrodynamics — insta- 
bility 



1. Introduction 



The mechanism of core-collapse supernovae has been an open question for decades. 
In the past 15 years, a lot of detailed numerical simulations have been done by many au- 



thors, paying attention to various aspects of the dynamics (see, for example, IfCotake et al. 



( 120061 ) for the latest review of the subject). Spherically symmetric models have incorpo- 
rated full general relativity and solved the r adiation-hydrodynamical equations with realis- 



tic equations of state and neutrino physics (jRampp &: Janka 1 12000| ; iThompson et al.l 12003 



Liebendorfer et ajj [2005| ; ISumiyoshi et al.ll2005l ). These computations demonstrated consis- 
tently that the stalled shock wave is highly unlikely to be revived by neutrino heatings in 
the spherically symmetric dynamics. 

Multi- dimensional s i mulati ons have been also making a remarkable progress over the 
times (see iKotake et al.l (120061 ) and references therein). Local non-spherical motion s of 



conve ction- type are common in the envelope as inferred from the spectra and light curves (INomoto et al 
20061 ). and are believed to occur also in the core. The global asymmetry as obs erved in the 



eject a of SN1987A is also supposed to be common to the core-collapse supernovae (jWang et al. 



20021 ). Motivated by these observations, the implications of the asymmetry for the explo- 
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sion m echanism have been investigated intensively and extensively (see again iKotake et al. 
( 120061 ) and references therein). 



Recently, the instability of the standing accretion shock wave was discovered by lBlondin et al 



(120031) in the context of core-collapse supernova and its nature has been studied both analyti 



cally ( Foglizzo et al. 



2006a 



Qhnishi et al 



2006a 



2006; Ya masaki fc Yamada 112006) and n umerically ( Blondin fc Mezzacappa 



Blondin fc Mezzacappa Il2006bl ). The ramifications to the global 



asym metry of the eject a as mentioned above and the high kick velocities of young pul 



sars ( ILyne &: Lorimer Ill994l ) were dis cussed by 



tions. Since it has been demonstrated (IHerant et al. 



Scheck et 



1994 



al. 



(120041) in more rea l istic simula- 



Janka fc Miiller 1996; Yamasaki & Yamada 



20061 ) over the years that these instabilities tend to fa cilitate the shock r evival by neutrino 
heatings, attempts have been done with little success (IBuras et al.l 120061 ) so far to show by 
realistic multi-dimensional simulations that an explosion is indeed induced. 

This situation may change. Very recently, Burrows et al. ( 2006a , b cl) have proposed a 
new possibility to induce explosions. They have done long-term simulations of axisymmetric 
2D radiation-hydrodynamics. Although some microphysics are incomplete and the formula- 
tion is essentially non-relativistic, their models with no equatorial symmetry and artificial 
inner boundary condition being imposed have obtained a shock revival after a rather long 
delay (> 500ms after bounce) not by the neutrino heatings but by what they call the "acous- 
tic" mechanism. They have claimed that the turbulence of accreted matter induces mainly 
g-mode oscillations in the proto neutron star, which in turn emit acoustic waves that are 
dissipated as they propagate outward to become shock waves. 

This is a quite new idea that has eluded the discovery by multi-dimensional numerical 
studies for many years, the reasons for which are, the authors think, that previous computa- 
tions were not long enough and/or put too restrictive conditions such as the artificial inner 
boundary t hat has been f r equent ly employed to reduce the CFL constraint on numerical 
time steps. iBurrows et al. (bo06aB 3) have insisted that the merit of this mechanism is the 
self-regulation in such a way that the g-mode excitations and subsequent emissions of acous- 
tic waves are maintained until the shock revival occurs. This is an interesting possibility and 
certainly worth further exploration. 

There are indeed a lot of things remaining to be clarified. First of all, the mechanism 
seems to be rather inefficient in transferring the energy. Inferred from the published re- 
sults, the explosion energy appears to be substantially smaller than the canonical value. 
Since the mechanism constitutes of multiple steps, the efficiency of the couplings between 
different steps should be understood. For that purpose, we think that the linear analysis 
utilizing an idealized model that mimics the realistic situation is more suitable than large 
scale simulations. In this paper, we calculate the eigen frequencies and eigen functions of the 
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g-modes in a non-rotating proto neutron star. Assuming that the standing accretion shock 
instability (SASI) is responsible for the excitation of the g-modes in the acoustic mechanism, 
we estimate the efficiency, formulating and solving the equations for forced oscillations of 
the g-modes. In so doing, we use the pressure fluct uations in SASI that we obtained in 
the previous numerical studies (jOhnishi et al.ll2006al ). The emissions of acoustic waves by 
these excit ed q- modes are curr ently being investigated numerically and will be published 
elsewhere ( jOhnishi et al.ll2006bl ). 

The organization of the paper is as follows. In section El we describe the models of 
proto neutron star and formulate the basic equations for the normal modes and their forced 
oscillations. The SASI that we use as an external force to drive the oscillations of g-modes 
is also presented there. We give the numerical results in section [31 We conclude the paper 
in section 0] with some discussions. 



2. Models &c Formulations 

We will investigate the g-mode excitations in the proto neutron star, one of the key 
steps in the acoustic mechanism, under the following assumptions: (i) A non-rotating proto 
neutron star is approximately in hydrostatic equilibrium, (ii) SASI in the accreting matter 
is responsible for the excitations of g-modes. (iii) The oscillation amplitudes in the proto 
neutron star are small compa red with the radius. As shown in the numerical simulations 
done by Burrows et al. ( 2006a b cl). these assumptions are reasonable after the stagnation of 



a prompt shock wave until the excited g-modes grow sufficiently and their non-linear effects 
on the dynamics become significant. Then the dynamics of the proto neutron star can be 
treated by the linear perturbation theory formulated below. 

The linear perturbation of stars can be described in terms of the Lagrangian displace- 
ment r), which obeys the equation of motion given as 

where C and F are, respectively, the self-adjoint operator and the external force driving the 
perturbations (see, e.g., Friedman & Schutz 1978; Unno et al. 1989). Here, we have assumed 
the unperturbed star to be static and spherically symmetric. Assuming the Lagrangian 
displacement to be 

t{t,r) = (2) 

where $, a and u a are the eigen function and eigen frequency, respectively, we can obtain the 
normal mode solutions, for which the eigen function and eigen frequency satisfy the following 
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equation, 

-a££ a (r) + C.* a (r) = 0. (3) 

Here, {a} denotes the mode index, which includes for spherical stars the angular eigen 
values £ and m (see below) and the radial mode number n. We impose for the normal mode 
solutions the regularity condition at the stellar center and Sp = at the stellar surface, where 
Sp stands for the Lagrangian change of the pressure p. 

When the unperturbed star is static and spherically symmetric, the orthogonality of the 
eigen functions allows us to adopt the normalization 

d 3 r p(r)CW ■ = MR 2 5 a , a , (4) 

where p, M and R denote the density, mass and radius of the unperturbed star, respectively 
(see, e.g., Unno et al. 1989). Any Lagrangian displacement £ that satisfies the same bound- 
ary conditions as those for the normal modes is then expressed as the linear superposition 
of the normal modes: 

ti(t,r) = J2Mt)Ur)- (5) 

a 

Substituting Eq. (jSJ) into Eq. p]) and using Eqs. and (J4]), we obtain the dynamical 
equation for the expansion coefficient A a (t) as 



d 2 A a (t) 
dt 2 



+ U ° Aa{t) = / rf3r PWCW-F(t,r). (6) 



2.1. Normal mode oscillations 

In this paper, we study adiabatic oscillations of a non-rotating spherical proto neutron 
star and employ the Cowling approximation, neglecting perturbations of the gravitational 
potential (see, e.g., Unno et al. 1989). For spherical stars, the separation of the angular vari- 
ables in the perturbation equations can be achieved by the use of the spherical harmonics in 
the spherical polar coordinates (r, 9, 0). The Lagrangian displacement vector £ is expressed 
in terms of the vector spherical harmonics as 

(£n to, = (W), rh(r)^, rh{r)-^ e ^j Y^A) , (7) 
and the Eulerian perturbation of pressure p' is given by 

p' = pgry 2 (r)Yr(e,<j ) )e^\ (8) 
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where g denotes the local gravitational acceleration. Eq. (j3J) is then reduced to 



^— = 0^ - 3) yi + [ - _ 2 -V g )y 2 



dz 



CiUT 



z Jl = ( Cl u 2 + tA) Vl + {-rA -U + l)y 2 , 
dz 



where 



V 9 

u 

Cl 
CD 2 



r/p, 

V 1 dlnp 

T r d In r 
dlnM. 



(ilnr 

{r/Rf/{M r /M) 

u 2 R 3 /{GM) , 
cHnp 1 dlnp 
dlnr Trflnr 



(9) 
(10) 

(11) 
(12) 

(13) 

(14) 
(15) 

(16) 



Here, T and M r stand for the adiabatic exponent and the stellar mass contained inside the 
radius r, respectively. The inner and outer boundary conditions are explicitly given by 



c\ijj 2 y\ — ly 2 — as r — > 



y\ - 3/2 

The total mode energy is defined as 



as r — > i? . 



(17) 
(18) 



1 
2 



y d 3 rp(r) ( 




2 

+ 




2 +(i: 


i 2 






dt 








Tp p 





(19) 



where p' is the Eulerian perturbation of density, and c s and N mean the sound velocity and 
the Brunt-Vaisala frequency defined, respectively, as 



c 2 = rp/p, N 2 = -gA 



(20) 



In order to obtain eigen frequencies and eigen functions, Eqs. (Q and (fTO]) are solved numer- 
ically as an eigen value problem for u), together with the boundary conditions (IT7|) and ffTS"]) . 
A Henyey-type relaxation method is employed (see, e.g., Unno et al. 1989). 



Introducing new variables £ and fj defined as 



£ = r yi exp 



zdr 



(21) 
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fj = gr 1/2 exp 



r N 2 
— dr 
o 9 



(22) 



we can rewrite the set of the pulsation equations and ( fTUj) in a canonical form given by 



(23) 



dry 



(^-iV 2 )f: 



dr r 2 H(r) 

where Lg is the Lamb frequency defined as 

Lg = c s 

r 

Here, the function H{r) is positive definite and is defined by 



(24) 



+ 1 



H(r) = exp 



N 2 g 



o V 9 



dr 



(25) 



(26) 



The canonical pulsation equations (1231) and (1241) are more appropriate to gain the insight 
into the qualitative properties of nonradial oscillations. Ignoring the spatial variations of all 
the coefficients in Eqs. f !23p and (I24p and employing a WKB-type local approximation, we 
obtain the eigen functions in the form, 



£(r), fj(r) oc exp(ik r r) 



(27) 



where k r means a radial wave number of the perturbations. Then the following dispersion 
relation is obtained: 



k 



10 c s 



(u 2 - Ll){uo 2 - N 2 



(2? 



which implies that if (u 2 — L 2 ){u 2 — N 2 ) > 0, the wave number k r is real and the corresponding 
wave can propagate in the radial direction, whereas the wave number k r becomes purely 
imaginary and the amplitude of the perturbation grows or damps exponentially with r when 
(u 2 - L 2 )(uj 2 - N 2 ) <0. 



2.2. Mode excitations 



In this study, we are concerned with the excitation of oscillations in the proto neutron 
star by the turbulence in accreting matter. In our formulation, we take into account the 
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latter effect as pressure perturbations at the surface of proto neutron star. Hence we set 
F = in Eq. (TjQ) and adopt the time- dependent non-vanishing pressure perturbation given 
by 

5p = f(t)p(r)Y e m (6,<f)) asr->R (29) 

instead of the ordinary condition 5p = at the proto neutron star surface. Here the function 
of time, fit), describes how the stellar surface is disturbed by the accreted matter and is taken 
from the numerical results obtained in our previous exploration of SASI (see section I2~3l) . 

Note that the boundary condition (I29p is no longer the same as that for the normal 
modes, 5p = at r = R. However, we can still use an expansion by the normal modes, 
£ a , for a part of the Lagrangian displacement by introducing a new dependent variable r\ 
defined as 



*7 = £ + C, (30) 

where 

(Cr,C*&) = (r (l) #1 (-^) V(M)> 0, O] . (31) 



Tj r=R £ + 3\RJ 



Substituting £ = r) — £ into Eq. ((T]), we obtain the master equation for r\ given by 

<9 2 T7 d 2 C 

8F + °-"=8F + c -<- (32) 



In Eq. (j3Tj) . C, is chosen so that rj could satisfy the same boundary conditions as those 
for the normal modes: 

5pr, = -— V ■ r] = at r = R , (33) 
P 

and rj is regular at r = 0. As a result, we can express 77 as a linear superposition of the eigen 
functions for the normal modes, 

r 1 (t,r) = J2Mt)Ur)- (34) 

a 

The expansion coefficient A a (t) is again determined by the equation similar to Eq. (J6]), 

^ + ^ -^pf^ CM ■ (§ + C ■ C) = ^jS. , (35) 
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where t = t(GM/R 3 ) 1 ^ 2 . In Eq. f l35l) . the driving-force term (R/GM 2 )S a is explicitly written 

as 



— Sa = am + b-^- 



(36) 



where a and 6 are constants determined by the background structure of the proto neutron 
star and the eigen functions for the normal modes and are given as 



--RJ0 



1 pR 3 



dx r ^ r x iJr2 
M 



Rp 



GMp \£ + 3 

2 ( 1 



Vg-l) V(y* - y{) 



x 



ci \£ + 3 



(rA-V g ) + l}yl 



--R JO 



M £ + 3 



Here the normalized radius x = r/R is introduced. 

Once we obtain A a (t), the Lagrangian displacement vector £ is given by 



(37) 
(38) 



(39) 



Since £ does not satisfy, in general, the orthogonality, 



Jd 3 rp(r)C(r)-t = 



(40) 



£ includes the components proportional to £ a , and the expansion coefficient A a (t) is only a 
part of the excited oscillation mode. The true amplitude B a (t) of the excited mode £ a in £ 
can be obtained by the inner product between £ and $, a as 



B a (t) 



' f <Pt p(r) &(t) ■ £ = A*® ~ I d 3 r p(r) C(r) ■ C ■ 



MR 2 J rw*«W * «W MR2 

Substituting Eqs. (EJ) and (!3~Pj) into Eq. ( |4B . we obtain 

B a (t) = A a (t)-bf(t). 



(41) 



(42) 



It is reasonable to assume that no modes are excited in the proto neutron star initially, i.e., 
B a (t) = and dB a (t)/dt — at t — 0. Accordingly, we start the integration of Eq. fl35l) 
with the following initial conditions, 



A Q (t) = 6/(t) , — ^ = at t = 



dt 



dt 



(43) 
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2.3. SASI 



As already mentioned above, we assume in this paper that oscillations in the proto 
neutron star are excited by the turbulence of accreting matter, which is supposed to be 
induced by SASI. SASI is a non-local hydrodynamical instability that occurs in the accre- 
tion flow through the standing shock wave onto the proto neutron star ( Blondin et al. 20031; 
Scheck et al.ll2004 ; Blondin fc Mezzacappa Il2006a ; Ohnishi et al.ll2006a i Blondin fc Mezzacappa 
2006bh . Although the mechanism is still controversial (jBlondin fc Mezzacappa Il2006al ). SASI 
is probably driven by the cycle of inward advectio ns of entropy- and vortex-perturbations and 



outw ard propagations of pressure fluctuations (jFoglizzo et al.l 120061 ; lYamasaki fc Yamada 
200d ). 



The prominent feature of SASI is the dominance of low I modes both in the linear and 
non-linear regimes. In particular, the I = 1 mode is thought to be a promising cause for 
the proper motions of young pulsars. In this paper, we regard SASI as a source of pressure 
perturbations exerted on the surface of the proto neutron star, that is, f(t) in Eq. (j2"9"]l . 



We employ the numerical results obtained by lOhnishi et al.l (j2006al ). in which we have 
done 2D axisymmetric simulations of the accretion flows throu gh the stand i ng sh ock wave 
onto the proto neutron star. The real istic equation of state by IShen et al.l (119981 ) and the 
standard reaction rates (IBruenn Ill985l ) for the neutrino absorptions and emissions on nucle- 
ons were employed. The mass accretion rate and the neutrino luminosity are the parameters 
that specify a particular steady accretion flow and were held constant during the compu- 
tations. The self-gravity of the accreting matter was ignored and only the gravitational 
attraction by the proto neutron star was taken into account. The mass of the proto neu- 
tron star was also set to be constant (1.4M ) so that we could obtain completely steady 
background flows, to which velocity perturbations were added initially. 

In the following, we use the results for the model with the mass accretion rate of 
M =O.19M /s and neutrino luminosity of L u = 3 x 10 52 erg/s, which are typical for the 
post bounce core around the shock stagnation. The neutrino spectra were assumed to 
take a Fermi-Dirac distribution with a vanishing chemical potential and a temperature of 
T Ve = 4MeV for the electron-type neutrino and T Ve = 5MeV for the electron-type anti- 
neutrino. The random velocity p erturbation of l ess tha n 1% was added to the radial velocity 
initially. We refer the readers to lOhnishi et al.l (j2006al ) for more detailed formulations. 



We observed SASI growing exponentially in the linear phase that lasts for ~100ms and 
then saturated to enter the non-linear phase. The I = 1 feature was clearly seen in the non- 
linear stage for this model although other modes were also generated by non-linear mode 
couplings. We expanded the Eulerian pressure perturbation at the inner boundary by the 
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spherical harmonic functions. We used only the I = 1 component in this study, since this 
was claimed to be the most important mode ( Burrows et al. 2006al lbllc1). 



The power spectrum of this component is shown in Fig. [TJ As can be seen, the pressure 
perturbation is actually a series of impulsive forces that last for quite short times. It should 
be noted that the linear analysis of SASI for this model gives the eigen frequency of ~ 30Hz 
for the linear I = 1 mode. The non-linear couplings broaden the spectrum both to the lower 
and higher frequencies. As seen later in section [21 the spectrum has most of the power in 
much lower frequencies than the typical g-mode frequencies of several hundreds Hz, which 
leads to a severe impedance mismatch in exciting (/-modes. 

It should be mentioned that the model studied here is not fully self-consistent. For 
example, what we need in Eq. f[2T?j) is the Lagrangian perturbation of pressure while what we 
obtained above is the Eulerian perturbation. This is simply because the former is difficult to 
obtain from the numerical data. Although they are different, rigorously speaking, we think 
that the difference will not change the outcome qualitatively, since it is the typical frequency 
and spectral feature of SASI that matters in this study. We should also bear in mind that 
the SASI spectrum employed in this study is just a single realization of essentially chaotic 
flows although we expect it is representative. 



2.4. Proto neutron star model 

Finally, the spherically symmetric, stationary proto neutron star that is adopted as the 
unperturbed state in this study is briefly described. The model is based e ssentially on the 



result s of detailed radiation-hydrodynamical simulation of core-collapse by ISumiyoshi et al. 



( 120051 ). but we have modified it in the following way. 



We have developed a numeric al code to compu te the long-term post-bounce evolution 



in the quasi-steady approximation (IWatanabe II2006I ). the details of which will be published 



elsewhere. In this formalism, we solve separately but consistently three different regions in 
the collapsing star simultaneously: (1) the proto neutron star that cools in a quasi-static way. 
We calculate both the equilibrium configuration and the neutrino transport in it with matter 
accretion taken into account. (2) the steady accretion flow from the standing shock wave 
down to the proto neutron star for the given neutrino luminosity obtained in the calculation 
(1) and matter accretion rate given by the next computation. (3) the dynamical implosion 
of the envelope toward the standing shock wave. The second and third regions are joined by 
the Rankine-Hugoniot condition at the shock wave consistently. 



The code can follow the spherically symmetric quasi-static evolutions for more than 10 



- 12 - 



seconds and we have confirmed that the evolution is well reproduced for the first second, 
for which the result of the detailed dynamical simulations is available. In this study, we 
utilize the proto neutron star, which is obtained in this approximation and is in hydro-static 
equilibrium at ~360ms after bounce. Physical quantities characterizing the proto neutron 
star model, that is, mass M, radius R, and Kepler frequency £\ are given in Table [TJ The 
profiles of density, pressure, specific entropy, and electron fraction are shown in Fig. [2] as a 
function of radius. It is noted that this spherically symmetric model contains a convectively 
unstable region around r ~ 20km. 

Here again, the proto neutron star model is not very consistent with the model for SASI. 
As mentioned in the previous section, it is not our aim to construct a fully self-consistent 
model in this paper, where the elementary feature of the g-mode excitations is explored 
qualitatively from an experimental point of view. For that purpose, the slight inconsistency 
in the current model does not matter. 



3. Results 

We first consider the oscillation properties of g-modes for the proto neutron star model 
described in the previous section. Before giving the eigen functions, we look into the so- 
called propagation diagram (see, e.g., Unno et al. 1989), which is useful to see where in 
the star the wave with a particular frequency u can propagate. The diagram shows the 
square of the Lamb frequency L 2 and the Brunt- Vaisala frequency N 2 as a function of the 
normalized radius r/R. As explained in section [2} since the waves can propagate when 
{uj 2 — L 2 ){oj 2 — N 2 ) > is satisfied, the region where this condition holds true is called the 
propagation zone, and the other regions in the diagram are referred to as an evanescent zone, 
in which the perturbations do not oscillate but grow or damp exponentially. 

Fig. [3] is the propagation diagram for the present proto neutron star model, in which 
we take i = 1 to calculate the Lamb frequency and the frequency is normalized according to 
Eq. (IT51) . The propagation regions can be classified into two categories: one is the g-mode 
propagation zone indicated as hatched regions and the other is the p-mode propagation zone 
shown as a cross-hatched region in Fig. [3j As mentioned already, since the proto neutron 
star has a convective region around r ~ 20km, which is an evanescent zone for g-modes, 
there are two disconnected g-mode propagation zones, one is near the center of the star and 
the other near the surface. This implies that a particular g-mode is almost trapped in one 
of these two disconnected g-mode propagation zones. In other words, the g-modes for the 
current proto neutron star model are essentially divided into two classes, i.e., what we refer 
to as the core g-modes (g c - modes) and the surface g-modes (g s -modes), whose amplitudes 
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of the eigen functions have a peak near the center and the surface of the proto neutron star, 
respectively. 

We now discuss the properties of the normal mode oscillations obtained numerically 
for the current proto neutron star model. In Table [2j we summarize the key quantities for 
several low-order g-modes: the normalized angular eigen frequencies, u, the corresponding 
frequencies in Hz, u, the mode energies, E, and the integral quantities, a and b, defined by 
Eqs. (157)) and (155)) . The last two are important in evaluating the excitations of g-modes 
by external perturbations in the next section. The mode energy gives us a rough measure 
of non-linear oscillations in the proto neutron star. Note that we employ the normalization 
condition for the eigen functions, Eq. (j3J), in evaluating the mode energy E. As expected, the 
mode energy is of the order of 10 53 ergs. The typical frequency of low-order g-modes is several 
hundreds Hz, a nd we can confirm tha t the numbers given in Table [2] are comparable with 



those found by iBurrows et al.l (l2006al ). It should be also emphasized that these frequencies 
are much larger than the typical frequency, <100 Hz, of SASI in the non-linear stage. It 
is also recognized in Table [2] that for £=1,2 the lowest mode is a surface mode and core 
modes and surface modes appear alternately as one goes to higher overtones. 

In Fig. HI the eigen functions, y\ and 2/2 in Eqs. (J7J and (j5J), are shown for the g-modes 
with i = 1 as a function of the normalized radius of the proto neutron star r/R. In the 
figure, the g n -mode means the n-th radial overtone g-mode, which is characterized by n 
radial nodes. As mentioned just above, the g±- and g 2 -modes are a typical surface and core 
g-mode, respectively. In fact, the eigen functions have a non-negligible amplitude either 
outside or inside the convection region, which is an evanescent zone for g-modes and located 
at r ~ 20km. This is also true for other g-modes with different I or n. Upon using these 
eigen functions, the quantities, a and b, which appear in Eq. (136)) and give the driving force 
for mode excitations, are obtained by integrating Eqs. (1571) and (1551) . Although the general 
trend is that their absolute values decrease as the radial node number increases, there are 
exceptions that are critically important in determining which mode is most efficiently excited 
by external perturbations as explained below. 

Now we move on to the excitation by accreting matter of the g-modes considered just 
above. As mentioned before, we assume in this study that the g-modes are excited by the 
pressure fluctuations around the surface of the proto neutron star, which are induced by 
SASI in the accretion flow. The pressure fluctuation f(t) in Eq. (1291) is obtained from the 
numerical results on SASI and is displayed in Fig. [1] as a power spectrum. It should be 
repeated that the major peaks of the spectrum are observed in the frequencies much lower 
than those of the low-order g-modes given in Table [2j In other words, the typical frequency 
of SASI, the driving force, is much lower than that of the low-order g-modes that should be 
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excited in the proto neutron star. It is true that these numbers are dependent on the models 
of proto neutron star and SASI as well as EOS employed, but we expect that the mismatch 
of the frequencies is a rather generic issue. 

By making use of the pressure fluctuations shown in Fig. [1] and solving Eq. (I35I) together 
with the initial condition given in Eq. (j4*3l) . we calculate the time evolution of the mode 
amplitude, B a (t), for several g-modes. Among them, the most efficiently excited g-mode is 
the g^-mode with I = 1. This is mainly because it has the largest a as given in Table [2J 
In Fig. 0, we show the temporal evolution of the mode energy for the g^-mode with i = 1. 
Since the driving forces, that is, the pressure fluctuations induced by SASI are a collection 
of impulsive forces, the mode energy is fluctuating on a short time scale. It is clear, however, 
that the mode is growing for ~100ms initially and then is saturated. The maximum value 
of the excited mode energy is < 10 50 (erg) for this mode, which is considerably smaller than 
the typical energy of non-linear oscillations, ~ 10 53 erg, given in Table [2j In this sense, 
the excitation of the g-modes is inefficient, at least for the present model. This is a direct 
consequence of the severe mismatch in the frequencies of the driving force, SASI, and the 
excited g-modes. Hence it is more appropriate to refer to what we see here as the forced 
oscillation rather than as the excitation. It is noted, however, that the Eu lerian pressure 



perturb ations at r ~ 50km are found to be of the same order as those found in lBurrows et al. 



f l2006af ). 

The maximum values of the mode energy for other g-modes, < 10 49 erg, are even smaller 
than that of the gf-mode with i = 1 for the proto neutron star model considered here. We 
show in Fig. [6] the temporal evolutions of the mode energy for the g\- and g2-modes with 
I = 1. The essential feature is the same as for the g2-mode shown in Fig. [51 These modes also 
grow for about 100ms initially and then get saturated with large fluctuations on a short time 
scale. It is emphasized again that all these excitations are non-resonant. In that case, the 
maximum amplitude of d 2 A a (t)/dt 2 can be approximated by the maximum value of af(t), 
which, for example, becomes ~ 10~ 2 for the g2-mode and an-order-of-magnitude smaller for 
other modes with £ — 1. From Eq. ffl9l . we then have approximate values of the mode 
energies, which are consistent with the numerical results given in Figs. [5] and [6j 



4. Summary and Discussions 



Inspired by a seri es of papers on the acoustic mechanism of core-collapse supernova by 
Burrows et al. ( 2006a b|,c]), we have calculated in this paper the eigen frequencies and eigen 
functions for several low-order g-modes in the proto neutron star, which were claimed to 
play an important role in the acoustic mechanism, and have studied also their excitations by 
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SASI in the accreting matter. We have constructed a proto neutron star model, which is sup- 
posed to represent a quasi-equilibrium configuration around a few hundre d milliseconds af- 



ter bo unce, based on the detailed radiation- hydrodynamical simulations by ISumiyoshi et al. 



( 120051 ). We have employed the temporal variation o f pressure on surface o f the proto neutron 



star, which were obtained in the study on SASI by lOhnishi et al.l ( ]2006al ). 



We have found that the g-modes are classified either as the surface g-mode or as the 
core g-mode according to the profile of the eigen function, since there is a convective zone in 
the proto neutron star, which is an evanescent zone for the g-modes. The eigen frequencies 
of the low-order g-modes are several hundred Hz, which are larger in general than the typical 
frequency of SASI, <100Hz, in the non-linear regime. 

Assuming that SASI is the driving force for the g-mode excitations, we have solved 
numerically the equations for the time evolution of the g-mode amplitude with the pressure 
perturbations at the surface of the proto neutron star taken into account. We have observed 
that the g-mode amplitudes grow for ~100ms initially and then get saturated. For the proto 
neutron star model considered in this paper, the gf-mode, that is, the second overtone core 
g-mode with i = 1 is found to have the largest mode energy at saturation. Even for this 
mode, the mode energy, < 10 50 erg, is much smaller than the typical energy, ~ 10 53 erg, 
for non-linear oscillations, which can be estimated by imposing the natural normalization 
condition to the eigen functions. This apparent inefficiency of the g-mode excitations is a 
direct consequence of the mismatch of the frequencies of the g-modes and SASI. 

It is noted, however, th at what is imp o rtant f or the acoustic mechanism may not be 



the mode energy. In fact, the lBurrows et al.l (l2006af ) mentioned that the g-mode oscillations 



acted like a transducer that transformed the gravitational energy of accreted matter to the 
energy of out-going sound waves. In the present formulation, unfortunately, it is in principle 
impossible to obtain the energy that will be transferred to the out-going acoustic waves. 
However, we may be able to infer it as follows. In the growth phase that lasts for ~ 100ms, 
the energy of accreting matter, or the work done by the external pressure perturbation in the 
present model, is transferred to the mode energy, as is apparent in Fig. [51 After saturation, 
on the other hand, the mode energy is almost unchanged. In this phase, then, it may be 
that the energy of the accreting matter is emitted just at the same rate as the accumulation 
rate. If this is really the case, we can estimate that the emitted energy, possibly as out- 
going acoustic waves, will be ~ 10 51 erg/s, which may not be a bad estimate, provided, in 
general, the explosion is considerably delayed (>500ms) in the acoustic mechanism. It is 
obvious, however, that more elaborate formulation of the problem, which would include the 
out-going w aves, is required. Num erical approaches from an experimental view point will be 



also useful (jOhnishi et al.ll2006bl ). 
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There are many other improvements remaining to be done. As repeated in the text, we 
have not attempted to make the model fully self-consistent in this paper, since the qualitative 
understanding of the g-mode excitations by SASI is the focus of the paper. However, the 
dependence of the results on the proto neutron star model, EOS, the progenitor model that 
determines the accretion rate, and so forth should be investigated in addition to making the 
model more self-consistent, all of which are currently under way. 
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Table 1. The mass, radius and Kepler frequency for the proto neutron star model. 



M/Mq 


R(km) 


^k(Hz) 


1.492 


55.6 


1074 



Table 2. Key quantities for g-modes in the proto neutron star. 



1 


Mode 




v (Hz) 


E(10 53 erg) 


a 


b 


1 


9i 


2.92 


499 


5.70 


3.50E-03 


1.72E-03 




9l 


1.43 


245 


1.30 


2.13E-02 


-2.99E-04 




9% 


1.03 


176 


1.97 


1.67E-03 


-1.54E-04 


2 


9l 


3.15 


538 


7.51 


5.91E-03 


9.97E-04 




9l 


2.24 


383 


3.15 


2.02E-03 


-1.01E-04 




93 


1.72 


294 


7.11 


1.02E-03 


2.38E-04 
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Fig. 1. — Power spectrum of the fractional pressure variation at r = 50km, given as a 
function of the frequency. Here, the fractional pressure variation is expanded in terms of 
spherical harmonics, and the fractional pressure variation associated with £ = 1 and m = 
is displayed. 
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Fig. 2. — Profiles of density p (upper left panel), pressure p (upper right panel), specific 
entropy s (lower left panel), and electron fraction Y e (lower right panel) for the proto neutron 
star model considered in this study. Here, ks denotes the Boltzmann constant. 
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Fig. 3. — Propagation diagram for the proto neutron star model given in Fig. [2, in which 
the square of the Lamb frequency L\ for i = 1 and the Brunt- Vaisala frequency N 2 are 
shown as a function of the normalized radius of the star r/R. The frequency is normalized 
according to Eq. (|T5|) . The g- and p-mode propagation zones are indicated as the hatched 
and cross-hatched regions, respectively. 
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Fig. 4. — Eigen functions y± and ?/2 as a function of the normalized radius r/R for the 
(/i-mode (left) and (fe-uwde (right) with i = 1 for the proto neutron star model given in 
Fig.d 




Fig. 5. — Time evolution of the mode energy for the ^2- m ode with t = 1 for the proto 
neutron star model given in Fig. [2j 
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Fig. 6. — Time evolution of the mode energy for the gi-mode (left) and (73-mode (right) with 
i = 1 for the proto neutron star model given in Fig. [2J 



